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Abstract 

We formulate a microscopic theory of the decay of a compound nucleus through fission which generalizes earlier 
microscopic approaches of fission dynamics performed in the framework of the adiabatic hypothesis. It is based on 
the constrained Hartree-Fock-Bogoliubov procedure and the Generator Coordinate Method, and requires an effective 
nucleon-nucleon interaction as the only input quantity. The basic assumption is that the slow evolution of the nuclear 
shape must be treated explicitely, whereas the rapidly time-dependent intrinsic excitations can be treated by statistical 
approximations. More precisely, we introduce a "reference density" pp which represents the slow evolution of the 
nuclear shape by a reduced density matrix and the state of intrinsic excitations by a canonical distribution at each 
given shape of the nucleus. The shape of the nuclear density distribution is described by parameters ("generator 
coordinates" q), not by "superabundant" degrees of freedom introduced in addition to the complete set of nucleonic 
degrees of freedom. 

We first derive a rigorous equation of motion for the reference density pp and, subsequently, simplify this equation on 
the basis of the "Markov approximation". The temperature T which appears in the canonical distribution is determined 
by the requirement that, at each time f, the reference density should correctly reproduce the mean excitation energy 
at given values of the shape parameters q. The resulting equation for the "local" temperature T(q, t) must be solved 
together with the equations of motion obtained for the reduced density matrix R{q\,q2\ t). 
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1. Introduction 



There are many systems in nature which are characterized by more than one time scale. Let us assume, for 
simplicity, that there is one time-scale for the slow, usually collective motion (described by variables or parameters 
q) and one much shorter time-scale for the rapid motion of a large number of individual particles (with the degrees 
of freedom x). Whereas the slow motion must be treated explicitly, the rapid dynamical evolution can often be 
approximated by a thermodynamical or, more generally, by some statistical "ansatz". There is a vast literature on the 
theoretical description of such systems. One refers to them as "transport theories". 

Their intense study started in the \9 th century with the observation of the "Brownian motion" [Q]], i.e. the motion 
of a pollen or dust particle submerged in a medium of the randomly moving molecules of a gas or of a liquid. As an 
average result of the randomly distributed collisions, the "Brownian particles" are exposed to a "friction force" —yV(t) 
reducing their average velocity v^f), and to a dissipative force measured by a parameter D. The dissipation parameter 

dT. - — a 

D is determined by the time rate — (r(f) - r(f)l of the mean square deviation of the trajectories r(t) from the average 
dt v ' 

trajectory r(f), the mean values being performed over an ensemble of trajectories. 

It was for the model of Brownian motion that Einstein derived the simple classical relation J^] : 

k B T 
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(M=mass of the Brownian particle, r=temperature, &B=Boltzmann constant). The generalized quantum-mechanical 
version of the Einstein relation ( 11. Il l is called "fluctuation-dissipation theorem" [3]. It will be discussed in Section[5] 
of our manuscript. 

The standard form of transport theories is based on a Hamiltonian H of the form fllit] 

H = H s {q) + H B (x) + H SB (x,q) (1.2) 

Hs depends on the slowly time-dependent variables q, H B , often referred to as "bath Hamiltonian", depends on the 
rapidly time-dependent variables x, and the Hamiltonian Hs B represents the coupling between both kinds of degrees 
of freedom. 

In the case of the atomic nucleus and numerous other systems, the complete set of independent degrees of freedom 
is given by the rapidly time-dependent variables of the constituent particles only. Here, the slow part of the time 
evolution is produced by a slow collective component of the particle motion. As an example, a nuclear fission or 
fusion-fission process proceeds through the slowly varying shape of the nuclear density distribution. 
If, nevertheless, one chooses to describe such a system by a Hamiltonian of type ( 11.2l i. the shape variables q are not 
independent of the particle variables x but related to them by constraints. Thus, for establishing a Hamiltonian of type 
( 11.21 . one needs a theory of collective motion as a prerequisite. Usually, one adopts a method proposed by Bohm and 
Pines [6]. It is based on the constraint that the collective variables q are to be equal to mean values of certain multipole 
operators q{x), which depend on the variables x of the individual particles : 

q = Tr(q(x)p) (1.3) 

Here, p is the density operator of the system or some approximation of it. Strictly speaking, fluctuations around these 
mean values exceed the scope of this approach. 

We, henceforth, consider in particular the nuclear fission process. One of the intriguing aspects of this process 
is the appearance of large fluctuations [7] of quantum mechanical as well as thermodynamical origin, as exemplified 
by the large widths of the mass- and kinetic energy-distributions observed even at small excitation energies of the 
fissioning system. 

In order to remedy this difficulty of the existing dynamical fission theories, we apply the Generator Coordinate 
Method (GCM) which has been introduced by Hill and Wheeler in their well-known paper of 1953 [8]. In this method, 
the wavefunction *¥(xi, . . ., x A \ t) of the nucleus is expanded in terms of a set of basis functions (f> a {x\, ■ ■ ■ , x A \ q) 
according to : 



x ¥(x u ...,x A ;t) = ^ j J dqf a (q,t)<p a (xi,. . .,x A ;q) (1.4) 

The wave-functions <p a (x\, ■ . ■ , Xa\ <?) depend on the nucleon variables X\, ...,Xa and on a limited number of "generator 
coordinates" (g.c.) q which serve to describe the possible shapes of the nucleus. The subscript "a" defines the quantum 
states of intrinsic excitation which may occur for each given shape q. The important aspect of the ansatz ( 11.41 ) is that 
the shape parameters q are summed over. Thus, they do not represent "over-abundant" variables. 
So far, in actual applications of the GCM to nuclear fission and to nuclear collective excitations, only th e groun d states 
<pa (xi, . . . ,xa\ q) of the Hartree-Fock-Bogoliubov (HFB) method have been superimposed in ( 11. 4t I9,[lft 111. These 
ground states are defined as the solutions of the variational problem 

= s(K(i)\h - Aq - nn$ - n P 2\K(q)> = (1-5) 

In Eq. dl .5b . the quantities H, q, N and Z represent the total Hamiltonian, a set of multipole operators, the neutron- and 
proton-number operators, respectively. The Lagrange parameters A, yu„ and fi p are defined by the constraints 

(<p a(> (q)\q\<Pa (q)) = q (1.6) 

<^( ? )|#I<M<7)> = N (1.7) 

(<Pa (q)\Z\<PaM)) = Z (1-8) 



where and Z are the neutron and proton numbers of the nucleus, and q is a set of generator coordinates. As usual in 
the constrained HFB method, the Hamiltonian H appearing in (11.5b must be built using an effective nucleon-nucleon 
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interaction adapted to the description of the intrinsic structure of the mean-field states 4> ao (q), such as the Skyrme 111 2 1 



or the Gogny lfl3ll effective interactions. In Refs. T(J Til, the Gogny force, which is known to provide a good 



description of both the average nuclear mean field and pairing correlations was used. 

Omitting intrinsic excitations <p a {q) + <Pa (q) in Eq. dl .4t implies that only low-lying vibrational excitations of the 
system are taken into account on the way to fission or fusion. Nevertheless, many aspects of the low-energy fission 
process like the coexistence of a fusion and a fission valley [9] and the basic features of the fragment kinetic energy 
distribution ifToh have been rather successfully obtained in spite of this restriction. 

In our present work, we want to go beyond these calculations in essentially two ways : 

(i) We include intrinsic excitations, i.e. our Hill-Wheeler basis derived from constrained HFB calculations in- 
cludes intrinsic excitations <f> a (q) up to a given upper energy limit. 

(ii) We represent the system by the density operator p(t) instead of describing it by a single state vector *F(£). 

As excited states (f> a (q), we use the eigenstates of the one-body part H[{q) in the Wick decomposition of the Hamilto- 
nian H'{q) (see Eq. (fTTOV ). 

As the part //' (20) + H' (02) of the Wick decomposition involving the creation (H' (20) ) or annihilation (ft'®®) of two 
quasi-particles vanishes as a result of the variational principe ( 11.51 ). the Hamiltonian H'(q) reads : 

H'(q) = H ,(p \q) + H' in \q) + H' {1) (q) (1.9) 

The Hamiltonian H'( \q) comprises all the terms with one creation and one annihilation operator, and H' {T) (q) is 
defined as the sum of all terms involving four quasi-particle creation or annihilation operators in normal order. Hence: 

H'M)Uq) = {H' (0) (q) + H' in \qj)Uq) = E\ a ( q )Uq) (UO) 

The "vacuum energy" at a given q is given by 



with the corresponding vacuum state <p ao (q) (see Eq.( 11.5l )). A presumably preferable choice of (f> a<s (q) would be 
provided by the solutions of the RPA at a given value of the g.c. q. Given the fact that it is difficult to deter- 
mine these multiphonon states for a heavy nucleus, we shall use the states defined in Eq. fll.lOb which represent 
multi-quasiparticle states in a nucleus with the mean multipole moments q. 
Two well-known practical difficulties are inherent to the GCM : 

(i) The Hill- Wheeler states <p a (q) are not orthogonal to each other. 

(ii) There may occur a number of linearly dependent or almost linearly dependent states among the basis states 
<p a (q)- This is revealed by the fact that the overlap matrix 

Na xai (quqi) := (<M<?i)l<M<?2)> (1.12) 
may have a number of vanishing or almost vanishing eigenvalues. 
We apply two remedies to cope with these difficulties: 

1) We work with discrete values of the generator coordinates (g.c.) instead of the continuous range of values implied 
in Eq.( ll.4b : Each given type of g.c. q varying from a lower limit inf(q) to an upper limit sup(q) is replaced by a finite 
number of discrete values with a spacing Aq: 

q — > inf(g) + nAq < sup(g) (1-13) 



where the positive integer numbers n range from to (sup(q) - inf(q))/ Aq. The step length Aq must be chosen in a 
physically reasonable way 

dq^AqY, = Yj (L14) 
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As we incorporate intrinsic excitations "a" for each set of the g.c. q, we surmise that the required accuracy of 
description can be achieved with not too small values of the step length Aq. 
2) We either use a biorthogonal system of basis states <p a {q), <p a (q) with 



(<p ai (qi)\<Pa 2 (q2)) = <W<W (1-15) 

or we replace the non-orthogonal Hill-Wheeler states <p a {q) by a system of orthonormal basis functions i(f a (q). 

Our paper is organized as follows: 
The just mentioned alternative of either using "adjoint" (="dual") basis states <p a {q) as "bra-states" in all the matrix 
elements or of introducing an adequately chosen orthonormal set of basis states ifr a (q) will be dealt with in Section[2l 
In Section[3] we shall introduce our basic approximation pp(t) for the statistical operator p(t) of the total system 
and derive a rigorously valid equation of motion for pp(t) using the formalism introduced by Nakajima and Zwanzig 
13. Q. This rigorously valid equation cannot be solved as it stands. An important and well-known step in transport 
theories is the application of the so-called "Markov approximation". Applying this approximation, we arrive at an 
integro-differential equation for pp(t) which is amenable to a numerical solution. 

The physical content of the Markov approximation is that the slowly time-dependent part pp(t) of the statistical 
operator p(t), which is designed to describe the slow collective part of the time evolution, should only depend in a 
negligible way on the detailed phase relations of the rapidly time-dependent intrinsic variables. The detailed form of 
the Markov approximation will be treated in Section[3] 

In the same Section, we shall also introduce a decomposition of the total Hamiltonian H into an average collective 
potential H co u and a remainder H cp i = H - H co u, and we shall postulate that the coupling Hamiltonian H cp i can be 
treated in perturbation theory. 

Within the slowly time-dependent part pp(t) of the density operator pit), the distribution of the system over the 
intrinsic excitations at given value of the g.c. q will be described by a canonical distribution which depends on a 
temperature T. This temperature T(q, ?), which should be considered as a useful parameter, will be chosen as a 
function of the generator coordinate q and of the time t in such a way that the ansatz pp(t) represents optimally the 
slowly time-dependent part of the density operator pit). The choice of the temperature T(q, t) and the equation by 
which it is to be determined will be presented in Sectiona l 

It is a well-known feature of the Generator Coordinate Method (GCM) that matrix-elements of the Hamiltonian 
or of other observables O between Hill- Wheeler states fll (q\) and (p^iqi) contain a factor which depends in a narrow, 
approximately Gaussian form on the difference r-q\ - qi, whereas the remaining part of the matrix element is a 
smooth function of the g.c. q\ and qi. Writing the elements of the overlap matrix in the form 

N aia2 (qi,q2) = {(f>aMl)\<Pa 2 (qz)) = e~ 1 n aiai (q\,q z ), (1-16) 

the factor n a ^{q\,qi) depends smoothly on the variables q\ and q2- Analogously, it is useful to define a "reduced 
matrix element" (<p ai {qi)\o\(f) a ,{q2)) by the equation 

, , rl , , „ (<Mffi)ld|<Mff2)) n 
(<M2i)|0|<M#2)> := — : (1.17) 

N ma2 (qi,q2) 

It is a convenient and usually good approximation [16] to assume that the reduced matrix-elements {<j) ai (q\)\o\4> ai {q2)) 
= (<p a ,(q\2 + f/2)\o\ (paXqn - r/2)) be represented by a Taylor expansion around q\ = qz = qn — iqi + qz)/2 up to 
second order. 

In transport theories, this approximation is customarily called the "Fokker-Planck approximation (FPA)". It was in- 
troduced in order to replace the "master equation" which is an integro-differential equation for a classical probability 
distribution depending on slowly varying variables, by simple differential equations of second order, the so-called 
"Fokker-Planck equations (FPE)". A pertinent presentation of this procedure is given in Ref. O- 

Section[5] of our paper will be devoted to the fluctuation-dissipation relation: It will be shown that, in the frame- 
work of our approach, such a relation holds for those transport coefficients in the NZ equation which are multiplied 
with the "diagonal" part R{q\,q\;t) of the reduced density matrix R(qi,q2', t) (see Eq. (15.101) ). Furthermore, as will 
become clear, the derivation of the fluctuation-dissipation relation relies on a few additional approximations beyond 
the mere perturbative treatment of the coupling Hamiltonian H cp \. 
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As we take into account the occurrence of intrinsic excitations of the fissioning or fusing system, our theory may be 
applied up to moderately high excitation energies. The upper limit of its applicability will be given by the appearance 
of final channels with more than two final composite particles, in particular, the "incomplete fusion" channels. In 
these reactions, only a part of the incoming projectile nucleus fuses and, subsequently, undergoes fission, whereas the 
remaining part is immediately leaving the interaction zone in predominantly forward direction. 
The emission of neutrons, protons, and photons by the fissioning nucleus will be described within the framework of 
our theory in a paper which we shall submit in due time to the same journal. For excitation energies in excess of 
about 50 MeV, the evaporation of neutrons and protons has been rather successfully described on the basis of the 
Thomas-Fermi approximation and a simple Fokker-Planck equation for the fissioning nucleus 11811 . 
Whereas the intrinsic excitations can only be treated quantum-mechanically, the collective motion leading to fission 
via the passage through strongly deformed nuclear shapes can approximately be described as a classical motion, 
at least at excitation energies where tunnelling is negligible. Indeed, valuable insights into the fusion- and fission- 
processes have been obtained within classical models 11911 . 

In our theory, the classical limit could be investigated by introducing the Wigner-Weyl representation for the statistical 
operator and by replacing commutators by Poisson brackets. We refrain from going into these details in the present 
paper. 

In Section [6] we discuss the relation of our theory to earlier work on the microscopic theory of nuclear fission. 
In Section|7] we summarize our results and point out open problems. 



2. Choice of basis functions 

The choice of the basis-functions is of great importance in our theory, because their dependence on the generator 
coordinates q enters our description of the slow collective motion and of the repartition of the intrinsic excitation 
energy.In the two subsections of this chapter, we deal with the biorthogonal basis and with a physically motivated 
special choice of an orthonormal basis. 

The Hill- Wheeler (HW) states introduced in Eq. (ll.lOt span a AA-dimensional vector space. If N gc is the number 
of discrete values of the generator coordinates and if one associates a given number Nmtr of intrinsic excitations with 
each given value of the generator coordinates, the dimension of the vector space is given by N — N gc • N m , r - 
The N Hill-Wheeler states <p a (q) and their adjoints <p a (q) defined by the property (11.151 l satisfy the completeness 
relation _ 

Yj \Uq))iUq)\ = Yj \Mq))(Uq)\ = i (2-D 

a,q a,q 

Henceforth, we convene to use the ordinary HW-states as "kef-states and the adjoints as "bra"-states. 

The adjoint states are obtained from the ordinary HW-states through the inverse (n~ 1 ) (qi, qi) of the overlap matrix 

(11.121 1, as can be easily seen : 

<M?T)I = 2< < M^)I'M^)X < M<72)I (2.2) 

a 2 q 2 

<^fe)I«M^)> = (AT 1 ) (quqi) (2.3) 

The labels a\, q\ characterizing a given adjoint state <p a ,(qi) denote the state with the largest amplitude (a^ 1 ) (q\,q2) 
in Eq. (12.21 i. We write the overlap matrix in the form 

N - I + K (2.4) 

The eigenvalues of the matrix K lie within the unit circle. Consequently, the series expansion for the inverse A^ 1 
converges 

AT 1 = (1 +K)' 1 = l-K + K 2 -K 3 + ... (2.5) 
The order of magnitude of the term K" in Eq. (12.5b is given by s" where the "overlap parameter" s is defined by 

s := exp(-A<? 2 /2o-o) (2.6) 



The parameter Aq is the distance between neighbouring discrete values of the generator coordinate q. To each type of 
generator coordinate corresponds a certain parameter s. 

In Fig. 1, we present the overlap matrix N for the simple case of 6 different HW-states. 



q& 



1 


N(q m ,q (2) ) 


N(q m ,q & ) 


N{q (2) ,q (l) ) 


1 


N(q (2 \q & ) 


N{q & ,q (l) ) 


N(q {3 \ q (1) ) 


1 



Figure 1: Form of the overlap matrix N for the case of 3 values (q 1 -^, q^ 2 \ q^) of the g.c. and 2 states of intrinsic motion (af \ a' 2 '). The 2x2 
submatrices for given values q^±q^ and q^—qU) of the g.c. are denoted by N(q('\ q^) and 1, respectively. 

The matrix-elements of a hermitian operator O - & between biorthogonal states satisfy the condition 

O aia2 (quq 2 ) := (<M&)\d\K(<l2)) = (<Pa 2 (q2)\0\Kiqx)r (2.7) 

The matrix O = {O aia2 (q\, q 2 )} is thus not hermitian. The lack of hermiticity of the biorthogonal matrix representation 
of a hermitian operator O is of no consequence for the thermodynamic mean value 

((d)) :=Tr (p(t)6) (2.8) 

given the fact that the trace is independent of the representation in which it is evaluated. 

Let us note that the adjoint basis-states <p a (q) have no physical meaning and that the operators 

fla(q) := \Uq))(Uq)\ (2-9) 

although satisfying the property 

n a (q)Tl a ,(q') = Tl a (q) (2.10) 

are not projection operators in the usual sense, as they are not hermitian. 

One can avoid the inconveniences of the representation in terms of biorthogonal basis functions by introducing a 
complete system of orthonormal basis states in the Af-dimensional HW-space. There are, of course, infinitely many 
possible choices of orthonormal basis functions. One possibility would be to use the orthonormal eigenstates of the 
overlap matrix N asa system of basis functions. 

A more physical choice is to determine the eigenstates of the hermitian operator 

a i q i a 2 qi 

where the hermitian matrix t H^ a2 (q\,qi) is defined by 

KaM'92) := \ (E[ ai {qi) + E' lai (q 2 j)N aia ,(q u q 2 ) (2.12) 

The energies E' la (q) pertain to states of independently moving quasi-particles in a nucleus whose shape is defined by 
the g.c. q (see Eq. dl.lOt ). The orthonormal eigenstates *¥ a of the Hamiltonian (12.1 j} with eigenvalues s a 

(^\H M \^ a ) = S a (^ a ) = BaSaf, (2. 13) 
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can be expanded in terms of the HW-states 4> a (q) or their adjoints 

Wa) = J] MArtXMlOWa) = J] \M))iM)\V a ) (2.14) 

a',q' a' ,q' 

In the limit of vanishing overlap parameter s, a given eigenstate TV becomes equal to a specific HW-state <p a {q). Thus, 
we can denote the orthonormal basis states *¥ a by T„(g) where the labels a, q characterize the main component <p a (q) 
in Eq. (12. 14-b which survives for e — > 0. 
The coefficients 

oZg := to&W (2-15) 
D%< := (frteW a ) (2.16) 



satisfy the relations 

YiBffi? = WfiWa) = S<# (2-17) 

aq 

^D^ q2 D^ iqi =S am 6 qiqi (2.18) 

a 

< ?1 = ((/> ai (qiWp) = Yj N ^(quqi)D; 2q2 (2.19) 

D°; qi = (^MxWp) = J] { N ~\ ia2 (Quqi)^ (2.20) 
Eq. (12. 13b can be written in the form 

2 0; i?1 *^C 2 («i-^)D? ?2 = bJS* (2.21) 

Ci?iC 2 92 

and the eigenenergies s a can also be obtained as solutions of the equation 

det{«£„(«i,«2) - e«iV aifl2 (?i,?2)} = (2.22) 

In the case that the overlap parameter s is noticeably smaller than 1, a given eigenstate Tv = T^g) of i/ M describes 
physically a state of independent quasi-particles whose main component is the eigenstate 4> a {q) of the HB-Hamiltonian 
H[(q). However, it contains admixtures of eigenstates <p a '(q') of the HB-Hamiltonian H[{q') with q' — q + Aq,q + 
2Aq, ... 

Let us note that the operator defined by 

tl a = |T„>(T ff | = Wa(q))(V a (q)\ (2-23) 
is a true projection operator satisfying not only the relations 

fijtp = 5 a pti a (2.24) 

but also 

tla = K (2.25) 

As most of the results to be obtained in the following chapters hold in the same form for the matrix-representation 
in terms of the biorthogonal or the orthonormal basis-functions, we introduce a unified notation encompassing the two 
alternatives in the following table 
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unified 
notation 


biorthogonal 
basis 


orthonormal 
basis 


\aq) 


\<p a (q)> 


!¥«(«)> 


(aq\ 






£' a (q) 




<(q) 


(aiqi\d\a 2 q2> 


<<M<7l)|0|<M<?2)> 





3. Derivation of an equation of motion for the slowly time-dependent density operator p>(f) and introduction 
of perturbation- and Markov- approximations 

The theory to be presented in this paper is devised to describe the slowly time-dependent average evolution of a 
fission- or fusion-process at moderate excitation energies. It should thus be applicable to reactions which proceed via 
the formation of compound nuclear resonances ("CNR") with a lifetime which is much larger than the time scale of 
its formation (< 1CT 22 s). The CNR are highly complicated nuclear many-body states which exhibit a finite width due 
to their coupling to the continua of the open decay channels. It is impossible to determine the detailed form of these 
states. According to the deep original insight of N. Bohr lEoll . the detailed structure of the CNRs is not important 
for the understanding of the slow time-evolution of the system in the fission channel which is the topic of the present 
paper. In the language of transport theory, Bohr's hypothesis implies that the "memory time" of the system is much 
smaller than the time-scale of the collective motion in the fission channel. Consequently, we may describe the slow 
collective time-evolution of the fission process starting from a simple initial condition for the collective flow. We shall 
come back to this question when discussing the "transport equation" we are going to derive. 

One usually considers the decay of an ensemble of CN resonances. But even in the case that we study the decay 
of a single, isolated CN resonance, the system should be described by the "density operator" (^"statistical operator") 
p(t) rather than by a wavefunction. The reason is that it is the density operator which expresses the information we 
have on a physical system and it is this quantity which simplifies when we ask for a reduced information only. 

The statistical operator satisfies the von Neumann-Liouville equation 

Pit) = j \H,p{t)] (3-1) 

In our case, H is the Hamiltonian of A nucleons (Z protons and N neutrons) interacting by effective nucleon-nucleon 
interactions v : 

6 =-t^ A ® + \t v ^ x » (3 - 2) 

1=1 i,;'=l 

The symbol x,- denotes the position, spin and isospin variables of the 2-th nucleon and M is the mass of the nucleon. 
We emphasize that the knowledge of the nucleon-nucleon interaction represents the only nuclear input information in 
our theory. 

The statistical operator pit) describes the entire time evolution, i.e. the slow collective change of the nucleonic 
distribution during fission as well as the rapidly changing state of internal motion which involves a superposition of 
many complicated intrinsic excitation modes. Of course, we also need a knowledge of the initial state p(0). Due to the 
effect of memory loss mentioned above, it is probably innocuous that we only have a rather incomplete information 
onp(O). 

At the origin of any transport theory, one has to introduce the part pp(t) of the total density operator which should 
describe the slowly proceeding dynamical evolution of the system only. As we do not have independent collective 
degrees of freedom at our disposal, we have to define the density operator pp(t) in terms of an "ansatz" for the matrix- 
elements {a\q\\pp{t)\a2q2) of pp{t) with respect to the basis-states introduced in Section 2. 

As the slow collective time change is not expected to depend on the precise form of the intrinsic motion, let us first 
define a "reduced density matrix" 

R{q u q 2 ;f)= p red {q x ,q 2 ;t) := ^(anlPiOWqi) (3.3) 

a 
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as the part of the operator pp(t) which is to describe the dynamics of the shape changes and, thereafter, of the motion 
of the fission fragments in space. 

Furthermore, we introduce a grand canonical distribution of the system with regard to the eigen-energies & a (q) of the 
basis states 

F a (q, T(q, t)) := — 1— exp (-/3(q, t)8' a (q)) (3.4) 

Zo(q, T) 



Here,/? is the reciprocal temperature 
and Z()(q, T) the partition function 



P(q,T)=-^— (3.5) 
T(q, t) 

Zo{q,T):=Yj e ~ PmG ' a{q) ( 3 - 6 > 



The temperature T may depend on the g.c. q and on the time t. The choice of the temperature T(q, t) as a function 
of q and t will be presented in Section 4. Let us only note that, for the case of a "nuclear transport theory", the 
"temperature" is not determined externally by the contact of the system with a bath, i.e. a large reservoir at a given 
temperature. Rather, we must consider it as a parameter which we have to choose so as to optimize the "ansatz" pp(t) 
for the slowly time-dependent part of the density operator p(t). 
We define the matrix representation of pp(t) by 

(aiqi\p P (i)\a 2 q 2 ) := S aia2 R(qt, q 2 ; t)F a2 (q 2 ,T(q 2 ,t)) (3.7) 

The Boltzmann distribution ( 13.41 ) which appears in the ansatz ( 13.71 ) is customarily also written in the form 

F a (q, T(q, t)) = exp [-ftq, t) (S' a (q) - G(q, T))] (3.8) 

where the Gibbs function G(q, T) is related to the partition function ( 13.6b by 

G(q,T) = -T(q,t)]nZo(q,T) (3.9) 

The following idea underlies the choice of the canonical distribution of the system with respect to the intrinsic excita- 
tions in Eq. ( 13.7b : 

As pp(t) should only describe the slow time-evolution of the system, the precise state of the intrinsic motion is not be- 
lieved to be relevant. The choice of the canonical distribution as approximation may be motivated by the observation 
that the distribution corresponds to the minimal information (^maximal entropy) on the state of internal motion (see 
Section 4). 

We draw attention to the fact that the matrix ( 13.71 ) does not define a hermitian operator pp(t), even if we use 
orthonormal basis states \aq). In addition, let us note that the reduced matrix R(q\,q 2 \t) is hermitian only if, in Eq. 
(13.3b . we use an orthonormal basis. It is not hermitian if we use the biorthogonal basis. A further discussion of the 
question of hermiticity will be given in the AppendixICl 

The non-hermitian character of pp{t) is of no consequence as far as mean-values of observables are con- 

cerned. We just have to evaluate the mean values with the hermitian part of pp{t) which is defined by 

p h P (i)= \(pp(t)+pl{tj) (3.10) 

We find 

((d)) :=Tr(p*(*)<5) = J] Re{<ai9i|0 P (f)|a2ft><fl2?2l<3|ai3i>} 

a l q l a 1 q 2 

=Re{Tr(pM0O)} (3.11) 
We also note that, for orthonormal basis states T^g), we can interpret 

R(q, q;t):= Yj^^WWaiq)) (3.12) 



as the probability to find the system in any one of the basis states f fl (g) for given q. In the case of the biorthogonal 
basis, the corresponding quantity 

R(q,q;t) := 2<W9)l0(Wa(?)> < 3 - 13 ) 

a 

cannot be interpreted as a probability, as \<f> a (q)){<f> a (q)\ is not a hermitian projection operator contrary to [¥ a (9))(^a(9% 
We can write the operator pp{t) as a projection operator P acting on pit): 

p P (t) = Pp(t) (3.14) 

Writing {aiq\\pp(t)\a2q2) as defined by Eq. ( 13.71 ) in the form 

{aiqi\Pp{t)\a 2 q 2 ) = J] (a 3 q3\p(t)\a 4 q4) (3.15) 

03^30494 

one finds that the matrix representation of the projection operator P has the form 

P ' a\q\aiqia^q^as,q\ — $ a\a%& qxq^& qiq$ a-sa$F a-SflZy T(q2-, 0) (3.16) 

Using the definition ( 13.16b of P we define a complementary projection operator 

Q = l-P (3.17) 

by the matrix representation 

Qa\q\a2q2^q^Aq4 — ^U: a? qi ^<XiM.\ ^qiq\ ^ asq-.ajqia-xq^a^qa. (3.18) 

One easily checks that the following relations hold: 

P 2 = P; Q 2 = Q; PQ = QP = (3.19) 

Through their dependence on the (slowly) time-dependent temperature, the projection operators are slowly time- 
dependent quantities. The total density operator p(t) can be decomposed into the slowly time-dependent part pp(t) 
and the remainder pQ(t): 

Pit) = (P+ Q)p(t) = pp(t) + p Q (t) (3.20) 

The part pg(f) of the statistical operator describes rapidly time-dependent processes and involves summations over a 
large number of terms with different phases. Whereas the detailed form of pg(f) will be considered in the next Section, 
let us note the properties 

J^(aq 1 \pQ(t)\aq 2 ) = (3.21) 

a 

Ti(p Q (t)d red ) = (3.22) 

Here, O red is an operator the matrix-representation of which is diagonal with respect to the quantum numbers of the 
intrinsic motion: 

(a 2 q2\0 red \a iqi ) := 6 a2ai O(q 2 ,qi) (3.23) 

Obviously, observables with the property (13.22b put to the test purely collective features of the system. Examples are 
the charge- and mass- numbers and the kinetic energy of relative motion of the fission fragments. 

For deriving an equation of motion for the operator pp(t), we make use of the formalism introduced by Nakajima 
1 1411 and Zwanzig B15I1 . Acting on the Eq. ( 13. U with the projection operators P and Q and taking into account their 
time-dependence, we obtain the coupled equations 

= - '-PL (p P (t) + p Q (t)) + P(t)p(t) (3.24) 

d -^- = - l -QL(p P {t)+p Q {t)) + &(t)p(t) (3.25) 
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where the Liouvillean operator L acting on an operator A is defined by 

LA := [H,A] (3.26) 

Using some properties of projection operators and their time derivatives, which are derived and presented in Appendix 
lAl we may write 

Pit)p(t) = P(t) (P{t) + Q(t))p(t) = P(t)p P (t) (3.27) 

and 

&{t)p(t) = -P(t)p(t) = -P(t)p P (t) (3.28) 
The coupled equations ( 13.24l >-( l3~2"5l i thus take the form 

d -PL=- l -PL (p P {t) + p Q (t)) + P(t)p P (t) (3.29) 
at n 

^-^- l -QLp Q {t) + &{t) (3.30) 
where the "source term" &(f) is defined by 

Ht) := - l -Q{t)Lp P {t) - P(t)pp(t) (3.31) 
h 

In order to obtain an equation of motion for the density operator pp(t), we have to find a formal solution of the Eq. 
(13.301 1 and substitute it into Eq. (13.291 1. Proceeding in analogy to Ref. 112111 . we introduce a Green function G(t, s) as a 
solution of the equation 

dGit, s) i s f „ 

— ^ + -Q(t)LG(t, s) = Q(t)6(t - s) (3.32) 
of h 

and write the operator pg(f) for t > in the form 



p Q (t) = G(t,0)p Q (0)+ f dsG(t,s)Q(s)fr(s) (3.33) 
Jo 

In the Appendix[A] it is shown that the operator P(t) has the property 

P{t)P(t) = (3.34) 

which implies 

P(t) = (P(t) + (2(f)) Pit) = Qit)Pit) (3.35) 
Consequently, the source term &is) satisfies 

Qis)fr(s) = &is) (3.36) 
and we may write, instead of Eq. (13.33b . more simply 



p fi (f) = G(f,0)p G (0) + f d s Git, s )fr(s) (3.37) 
Jo 

It is easily checked that the ansatz ( 13.371 ) together with the defining equation ( 13 .32b for the Green operator G(f, s) 
represents the solution of the Eq. ( 13 .30b with the initial value pg(0): 

dpQit) |337]<5G(f, 0) „ r° dGit, s) 



) r°° dGit, s) 

™ — — % <o, + r,s 

Jo 



-^Q(t)LG(t,0) + Q(t)5(t) 



- l -Qit)LGit, s) + Qit)6it - s) 
n 



Hs) 

dpoit) 13371 i - ~ - ~ «' - ~ r - i 

= - j.Qit)LGit, 0)p e (0) + e(0)p G (0)5(f) - ~<2(f)L [p e (0 - G(f, 0)p G (0)] + 6 it)Qit)&it) (3.38) 
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Apart from the term Q(Q)pQ(0)6(t), which expresses the initial condition for pg(f), the r.h.s. of Eq. (13 ,38b is seen to 
agree for t > with the r.h.s. of Eq. (13.30b . 

Substituting the formal solution ( 13 .37b into the Eq. ( 13.29b , we obtain the desired integro-differential equation for 
MO: 

= - l -P{t)W{t,Q)p Q {0) - l -P{t)Lp P {t) 

-~P(t)L I ds6(t,s)fr(s) + P(t)p P (t) (3.39) 
Jo 

We will refer to this equation as "Nakajima-Zwanzig (NZ)"-equation. It is a rigorously valid equation of motion for 
the slowly time-dependent part pp{t) of the total density operator p{t). 

The formal solution of the equation (13.32b for the Green operator G(t, s) can be written as the following time- 
ordered product 

G(t, s) = 6bC* " s)Q(t)t |«" -C dT ~ >(T) ]} . (3.40) 

T denoting the time-ordering operator. The detailed proof that ( 13.40b is a solution of Eq. ( 13.321 ) is given in the 
Appendix [B] Through its dependence on the Green operator G(t, s), the NZ-equation ( 13.391 ) still contains the full 
complexity of the total system. Nevertheless, this equation has the merit to display the different physical processes 
which contribute to the time evolution of pp{t). 

In order to describe the physical meaning of the different terms in the equation of motion ( 13. 39b . it is useful to 
decompose the matrix {a\q\\A\a 2 q 2 ) representing an arbitrary operator A into two parts using the projection operators 
P and Q: 

{aiq\\A\a 2 q 2 ) = (aiqi\PA\a 2 q 2 ) + (aiqi\QA\a 2 q 2 ) (3.41) 
<fli«i|M|0292> = 6 aia2 F 0t {q 2 ,T{q 2 ,i))A red (q u q 2 ) (3.42) 
A red (quq 2 ) := 2<fl?i£k»?2> (3.43) 

a 

(a iqi \QA\a 2 q 2 ) := {a x q x \A\a 2 q 2 ) - 8 a ^F ai (q 2 , T(q 2 , t))A red (q u q 2 ) (3.44) 

We shall refer to the two parts {a\q\\PA\a 2 q 2 ) and {ci\q\\QA\a 2 q 2 ) as the "canonical" or "P-part" and the "non- 
canonical" or "2-part" of the matrix {a\q\\A\a 2 q 2 ). 

Obviously, the Eq. (13.41b generalizes the decomposition of the density operator/) to one of an arbitrary operator A. 
Whereas {a\q\\PA\a 2 q 2 ) depends on intrinsic excitations only through a canonical distribution, the complementary 
part {a\q\\QA\a 2 q 2 ) is expected to depend on the quantum numbers a\ and a 2 in a complicated, possibly almost 
random way. 

The nature of the different terms in the NZ-equation ( 13.39b can be described as follows: 
Except for the term P(f)pp(f), all terms on the r.h.s. of the NZ-equation ( 13.39b contain the projection operator P(t) on 
their left, which means that their matrix representation is of the form (13.42b . (13.43b . The term P(t)pp(t) has the matrix 
representation 

(aiqi\P(f)p P (i)\a 2 q 2 ) = ^ Pa l q, ai q ia3 q 3 a i q i {ai ( li\Pp( t )\ a AqA) 

dF„,(q 2 , T(q 2 , tj) . 
=S aia2 g T T(q 2 , t)R(q u qi\ (3.45) 

where the derivative dF„ 2 /dT can be written in the form 

dF a2 (q 2 , T(q 2 , 0) F a2 {q 2 , T(q 2 , f)) f &' a Mi) ~ G(q 2 , T) 3G(q 2 , T)\ 

+ t-7. ' (3.46) 



dT T I T 8T 
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Introducing the entropy S(q2, T) of the system at the given value qi of the g.c. 



S(qi, T) = — (3.47) 

oT 

we can present the term {a\q\\P(t)p P (t)\a2q2) as follows 

(a iqi \P(t)p P (t)\a 2 q 2 ) =6 aiaz R(q u q 2 ; t)F a2 (q 2 , T(q 2 , 0)^|| ^ ~ ^ - S(q 2 ,T(q 2 ,t))\ (3.48) 

As the temperature is expected to be a slow function of time, the term is likely to be smaller than the other terms on 
the r.h.s of the NZ-equation. The role of the entropy will be discussed in Section|4] 

The inhomogeneous term -i/h P(f)LG(f, 0)pg(0) represents the memory of the system of the "non-canonical" part 
Pq(Q) of the initial density operator p(0). The operator G(t, 0) propagates p(0) from time to time t, so that 
-i/HLG(t, O)pg(O) essentially represents the time-derivative of G(t, 0)pg(0). Finally, the operator P(t) projects from it 
the component which is of P-type. One expects that the term -i/hP(t)LG(t, O)pg(O) fades away after a short "memory 
time" T mem < 10~ 22 s. In fact, the projection operator P filters from {a\q\\LG(t, 0)pq(0)|o2?2) the "reduced part" 

{a xqi \P(t)LG{t,Q)p Q (0)\a2qx) = S aia2 F a2 (q 2 , T) ^(aq^LGit, 0)p Q (0)\aq 2 ) (3.49) 

a 

As the density operator pg(0) represents the complicated part of the initial density p(0) and as the propagator G(t, 0) 
tends to complicate further this part of the statistical operator, the P-projection ( 13.491 ) is expected to become rapidly 
negligible. 

In all what follows, the inhomogeneous term -i/hP{t)LG(t, O)pg(O) of the NZ-equation will be neglected. The term 
was also neglected in all actual calculations of the fission dynamics we know of JD. H [l(| 181. 



Let us note that neglecting this term means essentially that one chooses the initial condition at a time where the 
collective motion towards fission has already started. The original compound system consists, at low excitation energy 
(< 8 MeV), of several narrow 1 eV) compound nuclear resonances and, at high excitation energy (> 10 MeV), of 
overlapping resonances of an average width of > 100 keV. Thus, just after its formation, the compound system does 
not at all resemble a canonical distribution over intrinsic excitations as required by pp(t), but rather a microcanonical 
ensemble of highly complex nuclear decaying states. The dynamical evolution of this initial state into one whose 
slowly time-dependent part is approximately canonically distributed has so far never been studied in detail. 

The terms -i/hP(t)Lpp(t) and -i/hP(t)L I ds ... in the NZ-equation d3.39t can be physically interpreted as 

Jo 

follows: 

As we can see from Eq. (13.29b , the term -i/hP(t)Lpp(t) would yield the whole time-derivative dpp(t)jdt if, at time 
f, the state of the system p(f) were entirely given by pp(t). However, there is also the part described by pg(f), and 
thus there is the contribution -i/hP(t)LpQ(t) to the time-rate dpp(t)/dt. As has been shown, this term takes the form 

Jr»oo 
dsG(t, s)d-(s). It describes that part of the rapid complicated time evolution which contributes to the 
o 

change of the slowly time-dependent density operator pp(t). As we shall see that this term describes friction and 
dissipation processes, we shall refer to it as the "dissipation term" of the NZ-equation. 

As it stands, the NZ-equation d3.39t cannot be solved due to the complexity of the propagator G(t, s). As this prop- 
agator appears in the equation of motion of the slowly time-dependent density operator pp(t), one may introduce some 
substantial simplifications, namely perturbation theory within the dissipative term and the Markov approximation. 

Let us first define a "collective potential Hamiltonian" H co u as a canonical average of the potential surfaces which 
correspond to the different intrinsic excitations of the system: 

In Eq. (I3.50l l, we first introduce a Hamiltonian H\ with the same eigenstates as the basis Hamiltonian HJ or resp., 
but with shifted eigenvalues 

Hi := ^ \aq) (8' a (q) + A(q)(aq\q\aq)) (aq\ (3.50) 



13 



or, separately for the biorthogonal and the orthonormal basis, resp., 

J] \cf> a (q)) (E"Z\q) + Mq){frM)\q\4> a (q))) (Mq)\ 
Hi - { * (3.50') 

= (&' a (q) + A(q)(aq\q\aq)) \aq) (3.51) 

The reason for shifting the eigenvalues by the amount A{q){aq\q\aq) is that H\ should be a simple part of the total 
Hamiltonian ft rather than of the Routhian ft' — ft - Aq. We denote the new eigen-energies by 

& a {q) := &' a (q) + A(q){aq\q\aq) (3.52) 

The matrix-representation of ft co u is then defined as 

(aiq\\H coU \a 2 q 2 ) = Y J ^a{q2)F a {q2,T{q 2 ,t)) (3.53) 

a 

i.e. as the canonical mean value of the Hamiltonian ft\. Due to the time-dependence of the temperature, the matrix 
{a\q\\H C0 [[\a.2q2) acquires a slow dependence on time. 

The eigen-energies & a (q) represent a family of potential surfaces, one for each quantum state a of the intrinsic 
motion. As we have discussed in Section 2, the quantum number "a" and the generator coordinate q design a particular 
HW-state <p a {q), if the biorthogonal basis is used, and they denote the largest component <p a {q) in the basis state ^aiq), 
if the orthonormal basis is used. 

In the limit of vanishing temperature, the canonical mean value of the different potential surfaces becomes equal to 
the lowest potential surface & ao (q)- 

Urn V & a (q)F a (q, T) = S ao (q) (3.54) 

T(q 2 ,t)->0 t—' 
a 

We define the coupling Hamiltonian H cp \ as the difference between the total Hamiltonian H and H c „u: 

H cp i = H - H co u (3.55) 

with the matrix-representation 

(aiqi\H C pi\a 2 q2) = (aiqi\H\a 2 q2) - {a\q\\H co u\a 2 q2) (3.56) 
Using the Wick decomposition for ft', we obtain a more explicit form of (a\qi \H cp i\a2q2)- 

(a l q l \H cp i\a 2 q2) =(aiqi\6 - A(q 2 )q\a 2 q2} + A{q 2 ){a^q\\q\a 2 q2) ~ {a\q\\H coU \a2q2) 

or 

{aiqi\H cp i\a 2 q2) =S aia2 S qiq2 (S' a2 (q 2 ) - ({B'iqi)))) + A(q 2 ) ((a l qi\q\a 2 q2) - ((<?(«z)») + («i?il-^2 (92)l fl 2?2> (3.57) 
Here, we have introduced a short notation for canonical mean values 

«£'(<?)» = J] K(<l)F a (q, T(q, f)) (3.58) 

a 



The first term in Eq. (13.571 1 represents the fluctuation of the intrinsic energies around their thermal mean value at a 
given value of the g.c. q 2 , and the second term the fluctuation of the "collective transport term" A(q2)(aiqi\q\a2q2) 

14 



around its thermal mean value at the g.c. qi. The last term in Eq. ( 13.57b is given by the matrix-element of the 2-body 
part ff 2 (#2) of the Hamiltonian H'(q2). The total matrix {a\q\\H cp i\a2qi) is hermitian if we use the orthonormal basis 
functions ^aiq): 

whereas the different terms on the r.h.s. of ( 13.571 ) are not hermitian. This is due to the fact that the Routhian H'(q) = 
H - A(q)q depends on the value of the generator coordinate. 

One can easily see that the Liouvillean L co u associated with H co u commutes with the projection operators: The matrix 
representation of the Liouvilleans read: 

^ai^ia2#2 a 3?3 a 4(?4 — ^2^4 ?2^4 

(aiq^Hfoqi) - S am 6 gm {a4q4\S\a 2 q2) (3.60) 

Lcplawawawatqi = ff a^«2?4< fl l«ll#*/l a 3«3> ~ Sa^Ms^mlHcpllaiqi) (3.61) 

{a A q4\H coU \a 2 q2) (3.62) 

(3.53) / \ 

«£(<?i)» - mm)))) 0.62') 

Using the matrix-representations ( 13.161 ). ( 13. 621 . one easily checks that the following relations hold: 

[P, L co „] = = [1 - Q, L coll ] = -[Q, L coll ] (3.63) 
Due to the relations ( 13.63b . the NZ-equation ( 13.39b can be written in the form: 

= ~ I P(t) (^'( f )G(f, 0)p e (0) + (L co „(t) + L cpl (t))p P (t)) - -^P(t)L cpl (t) j ds G(t, s)Q( S )L cp ,( S )p P (s) 

/-*oc 

+ l -P(t)L cpl (t) J o dsG{t,s)P(s)pp(s) + P(t)p P {t) (3.64) 

Eq. ( 13.64b is still rigorously valid. We shall now apply the "Markov approximation". 

This approximation is based on the observation that the non-negligible contributions to the integrals in Eq. (13.64b 
originate from a short time interval t - T re [ < s < t, where the "relaxation time" r re / is much smaller than the typical 
time scale r e of the collective motion. 

One may estimate T re i to be between 10~ 23 s and 10~ 22 s and r c to be of the order of 10~ 21 s, i.e. 

T re i « t c (3.65) 

The quantity pp{t) and the temperature T(q,t) (see Section [4} vary slowly, i.e. on the time scale t c , and so do the 
projection operators P(t), Q(t), and the Liouvillean L cp i(t), which depend on time through the temperature. 
The Markov approximation consists of several steps: 

1. The first step is to replace the time variable s by t in the slowly time-dependent quantities of the integrals: 

f dsG(t,s)Q( S )L cp ,( S )p P (s)^ f dsG(t,s)Q(t)L cp i(t)p P (t) (3.66) 
Jo Jo 

dsG(t,s)P(s)p P (s) I dsG(t,s)P(t)pp(t) (3.67) 
Jo 

2. Simplification of the Green propagator G(t, s) 

Replacing the slowly time-dependent quantities Q(f) and P(j) by Q(t) and P(t) in ( 13.40b . the time-ordering T 
becomes superfluous and we obtain 



--(?-.) 
G(t,s)^e (t-s)Q(t)e n 
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LQ(t) - -P(t) 
1 



(3.68) 



As the operator P(t) is proportional to f and as the temperature varies slowly as a function of time, the change 
of the temperature during the relaxation time, i.e. T re i ■ f(t) is a very small quantity. We, therefore, neglect the 
term proportional to P(t) in ( 13.68b altogether: 



- -t(t-s)LQ(t) 
G{t,s)^9 (t-s)Q{t)e h 



(3.69) 



The form ( 13.691 ) of the propagator is still too complicated for a practical application because it contains the 
Liouvillean L of the total system. 

There are two possibilities for further simplification: Either one considers the coupling Hamiltonian H cp i as a 
perturbation, or one approximates the propagator (13.69b by considering the matrix elements of H cp i as random 
numbers. The two ways of proceeding are in a certain sense complementary because the random matrix method 
implies that the interaction H cp i acts many times. 

In what follows, we shall use the perturbation theory for H cp i. In fact, perturbation theory underlies most of 
the work on transport processes and especially the work on nuclear fission Ji , 19]. The random matrix method 
was used by Weidenmiiller et al. [22] for formulating a transport theory of heavy ions reactions. It would be 
very interesting indeed to investigate the random matrix approximation of the propagator ( 13.69b parallely to 
the perturbation treatment. The formulation of the propagator ( 13.69b on the basis of the random matrix theory 
and the consequences thereof, for instance concerning the fluctuation-dissipation relation, would constitute a 
problem of its own merit, and exceed the scope of the present work. 

The perturbation theory implies that we replace the total Liouvillean L in Eq. ( 13.691 ) by the Liouvillean L\ 
corresponding to the Hamiltonian H\ introduced in Eq. ( 13.50b . One thus replaces the Green fuction ( 13.691 ) by 



G(t,s)^0 o (t-s)Q(t)e U 



(f-a)Lig(f) 



(3.70) 



As the basis functions \aq) are eigenfunctions of the Hamiltonian H\ (see Eq. (13.51b ), the propagator (13.70b 
becomes very simple if we may replace the projection operator Q(t) in the exponent of ( 13.70b by 1 : 



, --{t-s)U 
G{t,s)^e {t-s)Q(t)e h 



(3.71) 



The replacement of Eq. ( 13.70b by Eq. ( 13.71b would be strictly valid if the projection operators P{t) and Q{t) 
commuted with L\ . One easily checks the result 



(I ' ^-Dfli^i 



((y/ 2 (l;(/.!« 4 (/| 



(q2){p ai {q2) ■ 



& ai (qi) + & a ,(qi) - 8, 



; fl 3 (<72)) 



(3.72) 



It implies that the matrix element (fli^il [A^i] \aiqi), and consequently also (fli^il \aiqi) vanish for 

q\-qi- As the propagator G(t, s) appears in the "dissipative" integrals of Eq. ( 13.64b . the contributions of the 
matrix-elements (fli^i|G(f, s)|fl2<72) with q\-qi are expected to be much more important than the ones with 
q\+qi- As, within the subspace of matrix elements diagonal in the g.c, the Liouvillean L\ commutes with the 
projection operators, the approximation ( 13.71b of the expression ( 13.70b is acceptable. 

Introducing the approximations ( 13.66b , ( 13.67b . and ( 13.711 ) into Eq. ( 13.641 ). the NZ-equation assumes the form : 



at n 



-tt 



L cp i(t)e h p Q (0) + (L coU (t)+L cp i(t))pp(t) 



+ P(t)L cp i(t) 



P(t)p P (t) 



i r°° , --(t-s)L u 

J ds0 (t - s)Q(t)e ft L cp i(t)p P (t) 



i r°° , -Tit - s)L x x 

--J ds0 o {t-s)Q(t)e ft P(t)pp(t) 



(3.73) 



16 



Henceforth, we shall omit the inhomogeneous term in ( 13.731 1. Furthermore, in accordance with the approximation 
( I3.71l > of the Green function, we have put 

—z(t-s)Li . , -~(t-s)L\ 
Q(t)e h Q{t)^Q{t)e H (3.74) 

in the dissipative term of Eq. ( 13.731 1. 

In the last term of Eq. ( 13.731 1 we may leave away the projection operator Q(t) making use of the identity ( 13.35b . 
We can interpret this last term of (13.351 l as a modification of the dissipative term due to the time-dependence of the 
temperature. As the temperature changes slowly as a function of time, we expect this term to be small. 
We introduce the interaction representation of a Schrodinger operator (9(f) by 

Us-t)U^ l -{s-t)H x ^ - l -(s-t)ft x 
0\s -t;t)=efr 0(f) = eh 0(i)e h (3.75) 

The explicit time-dependence on t of the projection operators and of the Hamiltonian operators H co u and H cp i are 
produced by the variation in time of the temperature. Whereas the dependence of O'is - f; f) on the second argument t 
is slow, the dependence on (s—t) is rapid, because it is related to the intrinsic excitations of the system. It is satisfactory 
that the two different time scales are seen to enter the equation of motion in a natural and automatic way: 

d -^- = - l -P(t) [6 coU (f) + H cpl (t),p P (t)] + P(t)p P (t) 

i r°° 

- -^P(t) J ds 6 (t - s) [6 cpt (t), Q(t) [H' cpl (s -t;t), p' P (s - t; t)]] 

+ l -P(t) J ds e (t - s) [H cpl (t), e(f) e '' ( ' 5 -' )dl/fi >(f)pHf)e-' ( ' 5 -' )Al/;i ] 

Due to Eq. ( 13.34b and the relations P(t)Q(t) = and [Q(t), Li] — 0, we can simplify the two last terms of this equation. 
The NZ-equation for the density operator pp(f) thus reads: 



dpp(t) 
dt 



l -P(t) [H coU (t) + H cpl (t),pp(t)] + P(t)pp(t) 



i r°° 

- j^P(t) J ds 6 (t - s)H cpl (t)Q(t) [H' cpl (s - t; t),p' p ( S - t; t)} 

/-*oc 

+ -P(t) rfi6»o(?-*)^ ep /(?)e i( *"' )#l/ft ^WPp(Oe" i( '"' A/fi (3.76) 
n Jo 

It is straightforward to derive from Eq. ( 13.76b an equation of motion for the reduced density matrix R(q\ ,q2',t) defined 
by Eq. (13.3l i: Writing ( 13.76b in terms of a matrix representation and using (13.42b . we obtain 

8aia 2 F a2 (q2, T(q 2 , t)) ^(aqi\-^-\aq<i) =6 aia2 F a2 (q 2 , T(q 2 , t)) \ - l - ^\flq\\ \H coU (t) + H cp t(f),p P (t)] \aq 2 ) 

a \ a 

i r°° 

-tj ds6 (t-s)y^{aq 1 \H cp ,(t)Q(t)[Hl pl (s-tU),pU^- t ^)]\aq2) 
n Jo a 

+- J o ds6 (t-s) 2<^il^We ! ' (i " ,A/ ^(0pp(0e-'' (s - ,A/ft |a?2> J 
Summing over 01,02 and using the relation Y* a Fa(q,T(q,t))=l, we find the following equation of motion for the 
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reduced density matrix 

dR{q\,q 2 \t) i 



a'l ~R X {^313 .93)^(93, 92; t)F a3 (q 2 , T(q 2 , t)) - R(qi,q 3 ; t)F a3 (q 3 , T(q 3 , t))'Ha 3 aM^ 92)) 

a 3 q 3 

1 r°° ( 
- — I ds 6 (t -s) 2] ^cpi^iqm; m^^fe^^O?^; « - f, 0^(92, T(q 2 , r)) 

0»S?3?4 

- R 'a,aM^ S ~ f ' f )^ fl3 (?4, 7\<?4, t)) < H' cpla Jq 4 q2) 

+ ^ I dsG (t~s) e Hcpi aa (qiq3;t)F a (q 2 ,T(q2,t)) 

+ - dse {t-s)V9{ cp uAqiq3)K 3a Mm-,s-t,t)F a3 (q 2 ,T(q 2 ,t)) (3.77) 

Here, we have introduced the following notations: 

t Ka x a 1 {q\,q2) := (a l q l \H\a 2 q2) (3.78) 

^; aia2 (9i?2;0 := <ai?il^epz(0la2?2> (3.79) 

Klpi^qm) ■■= (a xqi \ff cpl {s - t; i)\a 2 q 2 ) = e /( ^.»2^'*)( fll?1 |H cp; (f)l«2^2> (3.80) 

R^iqm; s-t,t) := <fli ? i|-R 7 (i - f; r)|a 2 ? 2 > = <5 aia2 e i(i -^^ (9l ' ?2) /?(?i,?2; (3.81) 

R(t):= J] | fl y, >/?(?'„ ^;f)<a'?il (3-82) 

<a 1?1 |/?(f)|a 2 ft> = 6 aia2 R(q u q 2 ;t) (3.83) 
The frequency a> ai a 2 (q\,q2) occurring in Eq. ( 13.811 ) is related to the difference of eigen-energies (13.521 1 by 

fi^a t a 2 (qi,q2) := & ai (qi) - & ai (q2) (3.81') 
In the last term of Eq. ( 13.771 1. the factor F a% can be written in the form 

8F a Aq2, T(q2, f)) • 

F a M2, T(q 2 , f)) = 3 W ' W - t( g2 , t) (3.84) 

The equation of motion ( 13.771 ) for the reduced density matrix R{q\,q 2 \ f) must be solved together with the equation 
defining T(q, t) which will be derived in the next Section. 



4. Determination of the temperature 

So far, the temperature entered the "ansatz" for p pit) as an undetermined parameter. If we had exact solutions for 
pp(t) and pQ(t), the choice of the temperature would not matter, as it would only concern the decomposition of p(t) 
into pp(t) and /3g(f). Of course, in reality, we determine pp(t) from the equation ( 13.761 ). which was derived using the 
Markov approximation. One of the prerogatives of this approximation is that pp(t) should depend slowly on time. 
Even more so, we require physically that pp(t) should represent the slowly time-dependent part of the density operator 
pit) "as accurately as possible". It is this requirement which must determine the choice of the temperature. 

As the mean total energy of the system is independent of time, we should postulate in any case that the total energy 
of the system 

E:=Tr(6p(t))=Tr(6p h P (t)) (4.1) 
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is independent of time, where 

f> h P {t) = i{p P {t)+pl{tj) (4.2) 

is the hermitian part of p>. 

The condition (14. lb was used in Ref. J3l to determine the temperature as a function of time. In addition, at any given 
time f, the temperature T should also depend on the g.c. q, because the amount of intrinsic excitation is expected to 
vary along the dynamical paths of the system. In particular, we expect it to rise as the system approaches the scission 
region. Therefore, as a natural further condition to be fulfilled by the hermitian part p h p {t) of the reference density 
pp(t), we require that it should yield the same average energy as the total density pit) at any given value of the g.c. q. 
In order that this condition involves only real functions, it is given in the form 

1 2 Re {{aq\ [H,p(t)] + \aq)) = i J] Re (<*?! [ H ^"p (t )] + N>} (43) 

a a 

The symmetrized forms — [i/,p] + and — |i7,p^] + are used in (14.31 ) because the simple products Hp and Hp h p are not 
hermitian. In addition, we choose the real parts of the matrix elements because, in case of the biorthogonal basis, 
(<f> a (q)\ [#,p(r)] + \<p a (q)) and (<f> a (q)\ [#,p*(r)] + \(/> a (q)) are not necessarily real. 

Summing the Eq. ( 14.3b over the g.c. q and using the completeness relation ^ \aq){aq\-\, one finds that the Eq. ( 14.3b 

aq 

yields the condition ( 14. Il l of a conserved mean value of the total energy. 
The condition (14.31) can also be written in the form 



£ Re {{aq\ [H,p h Q (t)] + \aq)) = (4.4) 

a 

where Pg(f) is defined by 

p h Q (t):=p(t)-p h p(t) (4-5) 
For the sake of simplicity, in Eq. ( 14.4b , we replace the hermitian part Pg{t) by p@(f) : 

J] Re [(aq\ [H,p Q (t)] + \aq)} = (4.6) 

a 

If we use the orthonormal basis functions (\aq) = | v F a (q')), (aq\ = (^(g)!), the conditions ( 14.4b and ( 14.61 ) are exactly 
equivalent, as one can show that 

2 Re {<^(<7)| [6,p a Q (t)] + \V a (q))} = (4.7) 

a 

where 

P fl e (0 = ~(p fi «-p>)) (4.8) 

is the antihermitian part of pg(f). 

If we use the biorthogonal basis functions (\aq) = \(p a (q)), (aq\ = (0 a (<?)l)> tne condition ( 14.6b differs from ( 14.4b . 
because 

2 Re {(U4)\ [H,p a Q (t)] + \Uq))) * (4.9) 

a 

If the spacing Ag is such that the overlap parameter s <k 1, the expression ( 14.9b . which is of 0(b), is small. We thus 
expect that the simpler condition ( 14.6b yields almost the same temperature T(q, t) as condition (14.4b . Therefore, we, 
henceforth, use the relation (14.6b for determining the function T(q, t). 

Let us note that summing condition ( 14.6b over q leads to the conservation of the total energy 



Re (Tr(#p(?))} = Re {Tr(#p P (f))} = Re {Tr (Hp h p (t))} 
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(4.10) 



as we have 



Re {Tr (tip* (tj)} = ± ^Tt (Hp a P (t)) + Tr (#p£(*)) + } = (4.11) 
Let us now investigate the condition ( 14.6b in more detail: 

The equation (13.371) presents the density operator pgit) as a formal solution of the equation of motion ( 13.25b . Applying 
the Markov approximation ( 13.71b to the Green operator G(t, s), leaving away the memory term G(t, 0)pg(0), and 
replacing the source term &{s) in Eq. ( 13.37b by <x(f), we obtain the density operator /)(?(?) in the form 

i 

p e (f)- ds6 (t-s)Q(t)e h &(t) (4.12) 

Jo 

As for the condition ( 14.6b 

Re^^<a^| [#,p e (f)] + \aq) \ =Rej ^ ^{aq\6\aiqi){aiqi\p Q (t)\aq) 

+(aq\p Q (f)\aiqi){axq\\H\aq))] ) = (4.13) 

it is convenient to decompose pg(f) into two terms of different physical origin which correspond to the two terms of 
the source term (13.31b : 

Pe(0=P ( G 1) «+P ( e ) (0 (4.14) 

l :(t-s)L x 



PQ if) '- = ~nJ ds6 ^- s ^ e h 1 Q(t)L cpl (t)p P (t) (4.15) 

i 

PQ^^'nJ dsd ^- s ^ e h Lc P i(t)pp(t) (4.16) 

i 

X°° (f _ s )l. . 

ds0 o (t-s)e H P(t)pp(t) (4.17) 



Replacing the form (14.15b of p ( g(t) by the form ( 14.161 ) implies the same approximation which was used to replace the 
form ( 13.701 ) of the Green operator by the simpler form ( 13.71b . It is justified because, due to the factor {aq\H\a\q\), 
the main contribution to the condition ( 14.131 ) arises for q\-q. In the subspace of equal g.c, the projection operators 
commute with the Liouvillean L\, as one can see from Eq. ( 13.72b . 
After a straightforward calculation, one obtains 

• y-»O0 j 

(a\q\\p ( Q(t)\aq) =- -J ds8 Q (t-s) < {a x q x \ \H' cpl {s-t, t),p' p (s-t, f)] \aq) 

S aia F a (q, T(q, 0) Yj(bqi\ [H^is - f, t),p' P (s - t, ?)] \bq)\ (4.18) 
{a m \p { l\t)\aq) = -^ ds0 o (t - s)R' a Jq u q; s - t, t) 8 ^^' 0) f (g, t) (4.19) 

where H^^is - t, t) and p' p (s - t, t) are the operators H cp i and p> in the interaction picture (see Eq. (13.75b ). and where 
R' aia (qi , q; s - t , t) is defined by Eq. ( 13.81b . 

The temperature occurs in the condition (14.13b with different generator coordinates as arguments. As the matrices 
of p*Q and pfg are multiplied with matrix-elements 'Ha^iq, q\) and 'H aia (qi,q) which exhibit a narrow Gaussian-like 
dependence on q - q\, and as the "local" temperature T{q\,t) varies smoothly as a function of q\, we may introduce 
the approximation of replacing all the temperature functions appearing in Eq. (14.13b by T(q, t). 
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Let us introduce the matrices 7~ (1) and T~ (2) defined by putting all the T-functions occurring in the matrices p Q and 
pf by T(q,t) 

7<ll(quq; t) := [(a m \p%\aq)} w (4.20) 
Tf> a (quq;t)t(q,t) := {(am\pf\aq)} m (4.21) 

and correspondingly for the matrix-elements 7"^; (g, f) and T^](q, q\\ t). As can be seen from Eq. ( 14. 19b . the r.h.s. 
of (14.21b reads more explicitely 



T%l{quq;t) = {- f ds 6 (t - s)R ! aia (quq; s-t,t) 8Fa< ± T) ) (4.22) 

{JO 01 >T=T(q,t) 

The condition (14.6b thus assumes the form 

Re Yj r( \{f ?) + [W, ^ (2) ] + (f> ?) • *X«. 0) = (4.23) 

a 

or, more explicitely, 

Re | X (^.(^^^(Si.^O+^fe?!-.*)^?!,?)) 

+ 2 9", t) + Tg\{q,q l -,m ai a{quq))nq,f) \ = (4.24) 

The temperature condition (14.24b and the equation of motion (13.77b represent a coupled system of equations for the 
reduced density matrix R(q\,q2',t) and the temperature T(q,t). The solutions R(q\,q2',f) and T(q,f) are uniquely 
determined whenever initial values R{q\,q2\ t = 0) and r(^, f = 0) are defined. As an example, one could assume 
that the initial reference density p>(0) should correspond to a canonical distribution at the g.c. qo of the ground state 
valley. Then the initial values would be 

R{quqi\ 0) = S qm S q2qo (4.25) 
(aiqi\p P (0)\a 2 q2) = S qi (9o.r(? .0)) (4-26) 

The initial temperature could be defined by the knowledge of the total energy E m : 

^(ai^ol^|ai?o)^ a , (qo, T(q , 0)) = E tot (4.27) 
«i 

Let us comment the numerical problem of evaluating solutions of the coupled set of equations ( 13.77b for the reduced 
density R{q\,q2', t) and the equation (14.24b for the temperature T(q, t): 

Given the fact that a realistic description of the fission process requires a set of 2 to 3 generator coordinates, the 
technical problem is one of storage capacity or of rapid subroutines. As already mentioned, we replace the continuous 
values of the g.c.s by a network of discrete values and we also discretize the time. The step lengths Ag and At are to 
be determined by the required accuracy. 

From the knowledge of the initial functions R(q\,q2',t — 0) and T(q, t = 0), we first obtain R{q\,q2\t = At) from 
Eq. ( 13.77b and then T(q, At) from Eq. ( 14.24b . Subsequently, at the rfi 1 time step, one obtains R(q\,q2\ t„ + At) from 
substituting R{q\,q2\ t < t„), T{q, t < t„) on the r.h.s. of Eq. (13.771 ) and, analogously, one finds T(q, t„ + At) from Eq. 
(14.24b . If the storage capacity of the computer is sufficiently great, one can speed up the calculation considerably by 
storing the time-independent matrices {ci\q\\H\\a2q2), {a\qi\H C pi\a2q2), and the eigenvalues s a (q) before entering the 
routines for solving the Eqs. (13.77b and (14.24b . 

Through the Eqs. ( 13.77b and ( 14.24b the initial distribution R(qi,q2',t = 0) is thus propagated through the potential 
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landscape as a function of time up to the scission configurations which represent a (N - 1) surface in the space of N 
generator coordinates. The propagation of the reduced density R{q\,q 2 ; 0) up to values of R(q\,q 2 ; f) on the "scission 
surface" is reminiscent of a percolation problem. 

Let us conclude this Section with a few remarks concerning the entropy: 
As a first step, we have to introduce an adequate definition of the entropy. It seems to be appropriate to define a "local" 
entropy S {q, t) in relation to the "local" temperature T(q, t): 

S(q, T) := - Yj F a(<l> T ) lnF «(<l> T ) ( 4 - 28 ) 

S { q,T)= { ^ T))) T - G ^ T) (4.29) 

where 

«£'(<?, T))) := £>'a(q)Fa(q, T) (4.30) 

a 

The total entropy S(t) of the system is then to be defined by summing S(q, T) over all nuclear shapes weighted with 
the probability to find the system at the g.c. q at time f: 

S(t) r^Sfo.rfo.f)) *(«,$;*) (4.31) 

i 

We remind the reader of the fact (see Section 2) that R(q, q; t) can only be interpreted as a probability if we define it 
through the representation in terms of the orthonormal basis states *¥ a (q)- 

As long as we neglect the emission of light particles and photons during the fission process, we expect that the total 
entropy S (f) rises as a function of time 

dS(t) dS(q,T(q,t)) ■ 

_^ = ^ yq^yq, f) + ^ S (q, T(q, f)) R(q, qU)>0 (4.32) 

q q 

As the process of creating an increasing amount of intrinsic excitations, which are simply described by a canonical 
distribution, corresponds to a loss of information, the entropy must rise. However, as the functions T(q, t) and R(q, q; t) 
in ( 14.321 ) can only be obtained by solving the coupled equations ( 13.771 ) and ( 14.24b . we cannot prove the statement ( 14.321 ) 
in generality. 

In a realistic description of the fission process, deexcitation processes must be taken into account. As the emission 
processes reduce the entropy of the remaining system, the statement ( 14.321 ) can no longer be made. 



5. The fluctuation-dissipation relation 

The fluctuation-dissipation theorem relates part of the Fourier transforms of the response function and of the 
correlation function of a pair of operators with each other JUSIS]. The response and correlation functions appear in a 
natural way when studying the response of a system in thermal equilibrium with regard to small perturbations. As we 
chose to treat the coupling Hamiltonian H cp i in the integral term of the NZ-equation (13.771 ) as a perturbation, we may 
expect that some sort of fluctuation-dissipation relation should hold. As we do not deal with a thermal equilibrium, 
but with a slow passage of the system through thermal equilibrium states pertaining to different values of the g.c, we 
expect that a fluctuation-dissipation relation holds only in the limit of very slow collective transport. 
Apart from the last term in Eq. ( 13.77b . which depends on the time rate f of the temperature, the time integral in the 
Eq. ( 13.771 ) can be written as a sum of two terms 

£>(<?, , q 2 , t) := Y Diiquqi, q\, t)R(q[, qi\ + ^ D 2 (qi,q2, q\,q' 2 , t)R(q' v q' 2 \ t) (5.1) 
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where D\(qi,q2, q\,t) and £2(^1, #2, q\,q' 2 i are defined as follows 

I r°° 

D 1 (q u q 2 ,q[,i):=-—Y ds o (t - s)F a (q 2 ,T(q 2 ,t))e i< - s -' M ^ ) 

II aa >q' J ° 

X <?'; tYH^JLf, q[\s- t, t) - H cp i a , a ,(q U q'; tWc P iJq', q\\s- t, t)FAqi, T(q 2 , 0)} (5.2) 

i r°° 

D 2 (quq 2 ,q[,q 2 ,t):=-2Y ds 6 (t - s)FAq 2 ,T(q 2 ,t))e i ^-' )a} «' A ^ 

X frcpl aa ,{qi,q'vt) ( H I cpla , a {q' 2 , qi\ s-t,t)- ( Hc P i aa (quq[;t)'H I cpla , a ,(q 2 , qi\ s - t, t)F a (q 2 , T(q 2 , t))} (5.3) 

The matrix-elements of the coupling Hamimtonian H cp \ depend on the values of the generator coordinates in the bra 
and ket states through a narrow Gaussian (see Eqs. ( 11.161 ). ( 11.171 )). Thus, the matrix-elements have their largest values 
for equal values of the generator coordinates in bra and ket. Let us, therefore, consider the parts and D ( 2 0> of D\ 
and D 2 , which are of zeroth order in the overlap parameter s (see Eq. (12.61 )): 

i r°° 

Df\q u q 2 , q\,t) = Z \ ds 6 ° {t ~ S)F " iq2 ' T(q2 ' O^ - " 1 '^ 

x [Hcpi aa ,{qi,q\\tyH I cpla , a {q u qi\s - r,r) - i H cp i a , a ,(quqi;t)'H I cplaa (quqi; s - t, t)F a ,(q 2 , T(q 2 , r))} <W, (5-4) 
i r°° 

Df{ qi ,q 2 ,q\,q 2 ,t)=—Y ds 8 (t - s)FAq 2 ,T{q 2 ,t)y^- t ^'^ 

h aa> JO 

x [ H cpl aa ,{q\,q\\t) ( H I cp i a , a {q2, qi, s-t,t)- e Hc P l aa {qi,qiMYH 1 cpla , a ,{q 2 , q 2 \ s - t, i)F a (q 2 , T(q 2 , f))} <V,?i 8^ (5.5) 

The exponential factors in (15.2b to ( 15.51 ) (for their definition, see Eqs. ( 13. 81b . ( 13.811 ) refer to differences of energies at 
different q- values (qi,q 2 ) on the same potentiel energy surface characterized by given quantum numbers of intrinsic 
excitation "a" or "a"'. 

In the special case q\-q 2 , the exponential factors are equal to 1 and the coefficients assume the simple form 

i r°° 

Df\q u quq\,t) = J ds 6 (t - s){([9i cp i(qu qi) - (W f ((?i,?i)»„) 

■ {n'cpliVuqi) ~ «^(«i, 9 i)» gi )» e tf,^ (5.6) 

i r°° 

Df{q u q u q[,q' 2 ,t) = — J^ds 9 (t - s)(({<H' cpl {q u q x ) - ((Ti^foi, qi))) qi ) 

■ (U cp i{q U qi) ~ (( < Hc P l(quqi)))qi)))q 1 8q' l q l 8q' % q i (5.7) 

where Hcpi and < H 1 C ^ signify the matrices 

'Hcpiiquqi) := {ftcpia^iquqi-J)} (5.8) 
Hpfauqi) := \Hl plaia {q x ,q 2 - s - f,f)} (5.9) 

and where we denote the thermal mean-values of an operator O at a given value of q of the generator coordinates by 

{{0)) q -^YjFai^nqJMaqfilaq) (5.10) 

a 

The integrands on the right hand sides of the Eqs. ( 15.6b , (15.7b are seen to represent the correlation between the 
fluctuations 6H cp i(t), SH' cpl (s - t, t) of the operators H cp i(t), H ! cpl (s - t, t) around their respective thermal mean values. 
Decomposing the product of the fluctuations into a commutator and an anticommutator, we may write the coefficients 
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( 15.6b and ( 15.7b in the form 

i r°° 

Df(q uqu q\,t) = -— J ^e (f - «) {({[m P i(t), my,s - », *)]_», + - t, o] + » ? } ^ 

(5.11) 

Df\ququq[,q' 2 , t) = ^Jrf* Qo(* - s) - t, f), W cp /(0]_», + «[<W^(* - f, f), <V,?.<W 

(5.12) 

The mean-values of the commutator parts of Eqs. ( 15.1 U and (15.12b are proportional to the "response function" defined 
by — 

Xiq, s'; t) := l - ^ 6o(-s')F a (q, T(q, t)){aq\ [6H' cpl (s', t), 6H cpl (t)\ \aq) 



6o(-s')(([d<H I cpl ( S ', t), ffHcpM]_))q (5-13) 



h 

Here, we used the following definitions 

s' = s-t (5.14) 
6H' cpl (s', t) := H' cpl (s', f) - ({H ! cpl (s\ t))) q (5.15) 
6H cpl (t) := tf cp/ (0 - {{H cp ,(t))) q (5.16) 

The slow time-dependence on t is produced by the dependence of the temperature T(q, t) on t. The time-dependence 
on sf is rapid due to the dependence on the intrinsic excitations of the system. The coefficients D\ \q\, q\, q' v t) and 
D < ?{q\,qi,q' v q' 2 , t) being of very similar structure we, henceforth, shall refer to D^\q\,q\,q' v t) only, as the results 
obtained for Dj 0) can be immediately transferred to DfK 

The purpose of the following investigation is to show that a fluctuation-dissipation relation holds for these parts 
of the coefficients D\ and D^. 

To start with, we remind the reader of the fact that the coefficients Dj and are defined by the parts of order 
s° in the overlap parameter s. Consequently, by consistency, we have to put 

(aq\6H' cpl 5H cp i\aq) = ^(aq\SH' cpl \a'q')(a'q'\6H cp i\aq) 

a' V 

- Y J (aq\S6 I cpl \a'q}(a'q\86 cp i\aq) (5.17) 

a' 

Let us introduce the Fourier transform x(q, to; t) of the response function ^(^, s f \ t) 

X+oo 
ds' X (q,s';t)e^, (5.18) 
oo 

the so-called "dynamical susceptibility", and its "dissipative" partp^g, w; t) defined by 

X d (q, o;t) = - (xiq, or, t) ~x(q, -o)\ 0) (5.19) 
2i 

It is also given by 

X+oo 
ds'xiq, /;f)sin(ws') (5.20) 

It is easily seen thatp^ can be written in the form 

X d {q, co;t) = —J ds' J] F a (q, T(q, t))(aq\ [SH' cpl (s' , f), 6H cpl (t)] \aq)e^' (5.21) 
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The important point is that no Heaviside functions appear any more in the representation ( 15. 21b of the dissipative 
susceptibility. Due to this feature, the r.h.s. of ( 15. 21b is proportional to the Fourier transform of the antisymmetrical 
part A(q, s'; t) of the function 

C(q, s'; t) := J] F a (q, T(q, t))(aq\SH' cpl (s' , t)SH cpl (i)\aq) (5.22) 

a 

which represents the correlation between the fluctuations 8H cp i(t)=6H' cpI (0, t) and 5H I cpl {s' , f). 

Analogously to the decomposition of the expressions ( 15.1 lb . (15.121 ). we write d5.22l i in the form 

C(q, s'; t) = S(q, s'; t) + iA(q, s'; f), (5.23) 

where the symmetrical and the antisymmetrical parts of the correlation function C are defined by 

S {q, s';t):= ^ J] F a (q, T(q, t)){aq\ [6H' cpl (s', t), 6H cpl (t)] + \aq) (5.24) 

a 

A(q, /; f) := i ^ F a (q, T(q, t)){aq\ [sff cpl (s', f), 6H cpl (t)] \aq) (5.25) 

a 

The Fourier transforms, sometimes referred to as "spectral functions", are given by 

ds'C(q, s'; t)e i0JS ' = S(q, w; f) + iA(q, u; t), (5.26) 

oo 

where S (q, cj; t) and A(q, a>; t) are the Fourier transforms of S (q, a>; t) and A(q, a>; f). Comparing the r.h.s. of Eq. (15.21b 
with the definition of A(q, to; t), one obtains the relation 

X d {q,^t)= l -A{q,a)-,t) (5.27) 

which is one way of writing the fluctuation-dissipation relation. 

In what follows, we shall show that the r.h.s. of Eq. ( 15.27b can be written in the form 



where f3(q, t) is the reciprocal temperature 



l -A(q,cj;t) = ^-(l-e-^)c(q,cj;t) 



(5.28) 



P(q,t):=— 1— (5.29) 
T(q, t) 

The proof of the relation ( 15.28b is based on the "stationarity" of the mean-value 

((6ff cpl (s', t)5H cpl (t))) q = ((6H cpl (t)6H ! cpl (-s\ t))) q (5.30) 

and on the following symmetry relation for the correlation function 

C(q, s'; f) = C(q, -s'- ih/3\ t) (5.31) 

Although these two relations are well-known (see for instance Ref. fl), we have to discuss their validity, which is not 
automatic in the context of our theory: 

As a prerequisite of the proof, we have to accept the approximation ( 15.17b of only considering terms of order e°. We 
can thus use the cyclic invariance of the trace in a first step 

((6H' cpl {/, t)6H cpl {t))) q :=2><?l 7 , ; e * S Sl SH cpl (t)e-* s H '5H cpl {t)\aq) 



Zo(q) 

Z, - e -/%)#;(?) . , 

(aq\e-* s H, 6H cp ,(t) e* s 6H cp ,(t)\aq) (5.32) 

A>{q) 
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ans, subsequently, the commutation property 

[6u6[(q)]=0 (5.33) 
in a second step, for rewriting the r.h.s. of Eq. ( 15.321 ) 

_ -/%)#;(?) 
= Y^aqW^i-s' ,t) — SH cpl (t)\aq) = ((6H cpl (t)5H' cpl (-s' ,t))) q 

This proves relation ( 15.30b . 

Let us now consider the explicit form of the function C(q, —s'— iftfi; t) 

{aq\ - - - S qVtl SHU-s', t)e- MH ^H cp ,(t)\aq) (5.34) 

a z o(<?) 

If the Hamiltonian H' { (q), which appears in our equilibrium distribution, were equal to the Hamiltonian H\ defining 
the interaction picture, the r.h.s. of Eq. ( 15.34b would be equal to C(q, s', t) because of the cyclic invariance. 
In our case, these two Hamiltonians differ because we do not consider a true thermal equilibrium but a constrained 
one. In the expression ( 15.341 1. we may replace the operator H\ (see Eq. d3.50'l l) by H\{q), i.e. by the part of H\ acting 
on the subspace of states \aq) with fixed value q : 

H\(q) = 2j \ a< l) (fifl(?) + Mq){aq\q\aq)) {aq\. (5.35) 

a 

The Hamiltonian H[ (q) 



H[(q) = J]\aq)S' a (q)(aq\ (5.36) 



differs from Hi (q) by A(q)q(q) 
q(q) being defined by 



Hdq) = H[(q) + A(q)q(q) (5.37) 
q{q) ■= ^ \a'q)(a'q\q\a'q)(a'q\ (5.38) 

a' 

Eq. (15.34b can thus be written in the form 

C(q, -s'- ihJ3; t) = Y{aq\-^e mA(q) « (q) 5H'(-s', t)e^^^ lt ^^SA c ^\aq) 

a Zo(q) 

(aq\6H I c J-s\t)—-—e-P (q)A( ^ (5.39) 

a z o(q) 

In Eq. ( 15.39b , the matrix-elements 

(a'q\e- mA(q)q(q) 6H cp i(t)^ 

appear. The symmetry (15.31b of the correlation function C then holds, if we may put the exponential factor approxi- 
mately equal to 1 

e -P(q)Kq)((<*' q\q\a' q)-(aq\q\aq)) ^ y « ^q-j 

The approximation ( 15.40b holds, if the multipole moments in different intrinsic excitations do not fluctuate too much. 

Let us assume that the approximation ( 15.40b is valid. We then claim that the following relation holds for the 
Fourier transforms of the correlation function 

C(q, to; t) = C(q, -a>; t)S q)n *> (5.41) 
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as a result of the symmetry property ( 15.3 It . The proof is straightforward: 



- f + °° ■ - I53T1 f +0 ° 

C(q,to;i) = ds'C(q,s';t)e ,0JS U ^- J ds'C(q,-s'- ihj3;t)e" 

•J OO %J —CO 



L 



ds"C(q, s"; t ) e - ibJS "e mh ^ = C(q, -co; t)e mho 



As an immediate consequence of the relation (15. 41b we obtain 

S(q, to; t) = \- (l + e-#<? )fe ") C(^, -w; f) (5.42) 
A(q, to;i)=—[\- e - mhoJ ) C(q, -to; t) (5.43) 



2/ 

Substituting the result ( 15.431 1 into Eq. ( 15.271 1. we find the conventional form 

X d {q, to;t) = ± h (\- e~») C(q, -co; t) (5.44) 

of the fluctuation-dissipation relation 01. 

In Section Q] we had mentioned the simple classical result of Einstein (see Eq. ( 11. U ) relating the diffusion coef- 
ficient D to the friction constant y for the case of Brownian particles of mass M immersed into a gaseous or liquid 
medium. 

The classical limit of the relation (15.44b is obtained in the limit 

p{q,t)tko=^-^\ (5.45) 
T(q, t) 

In the approximation 

e -pf*> ~ i -pfta (5.4 6 ) 
the fluctuation-dissipation relation (15.44b assumes the form 

C(q, -co; t) * — x d (q, to; t) (5.47) 

CO 

It does no longer contain explicitly Planck's constant and corresponds to the Einstein relation ( 11. U . 
Using Eq. (15.42b . one can also relate the Fourier transform S (q, u>; t) to ~x d {q, co; t): 

S(q,co;t) = hcoth ^ q, ^ h0) y?{q, to; t) (5.48) 

with the classical limit 

S(q,co;t)^ — ^(q, to;f) (5.49) 
co 

It is thus shown that the parts of order O(e ) of the coefficients D\ and D2 (eqs. (15.2b . (15.3b ) do fulfill a fluctuation- 
dissipation relation. This is expected from the fact that in the limit O(e ), the coefficients D\ and D2 do not describe a 
transport of the system between different values of the generator coordinates but rather a quasi-stationary equilibrium 
at a given value of the generator coordinates. 

6. Simple limits of the theory 

Let us consider the equation of motion ( 13.77b for the reduced density matrix R(qi,q2; t) in some special cases: 
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If we restrict the Hill- Wheeler space to the grounstates <p ao {q) at each given value of the g.c. q, eq. (13.77b assumes 
the simple form 



dR(q\,q2',t) i 
dt h 



E [%oao(?l> 93)^(93. «25 ?3! O^Jdoofe, 92)] (6-1) 



This equation is equivalent to the description of the fission process by the usual GCM without intrinsic excitations. 
The latter theory underlies the calculations of Refs. ifl fioll . Using this theory, it was found that, at small excitation 
energies (< 10 MeV), the calculated mass distribution of the fragments in 238 U agreed quite well with the experimental 
one, although its width was a little too small ifioll . In this respect, it should be noted that the omission of intrinsic 
excitations does not mean that vibrations of the nucleus shape on the way to scission are neglected. On the contrary, 
the dynamical role of these vibrations was proved to be important to explain the large width of the calculated mass 
distribution. 

Another approximate version of the equation of motion ( 13.77b is obtained in the case that the "local" temperature 
T(q, f) is small compared to the energy of the lowest elementary excitation a + a^. 

T(q, t) « Min (S' a (q) - S' ao (q)) (6.2) 
The Boltzmann factors F a (q, T) can be written in the form 

e S'MIT AVaW-V^IT 

F a (q, T) = = ; — (6.3) 

Za(q,T) l+X^ ao e< E W- s ''iMl T 

which, in the limit of Eq. (16.21 i. yields 

F a (q, T) - 6 aao (6.4) 

Using the same technique, it is easy to verify that the time-derivative F a -t dF/dT goes to zero in the same limit. If 
we substitute these limiting values into (13.771 1. we obtain 

~~~7T^ = ~ J. E fa<° a ° (9i- 93)^(93, qi\ t) - R{q\,qi\ OH^ofe, qi)\ 

- ¥ ds e Q (t - E E {Hcniajmqr, W^iqw, s - 1, o^W^x-o/?^ 92; 

93 94 a+aa 

~ Mcpl^iqW, m l cplaaa {q Aqi ; s - t, f )e*W**X*-«) R{q ^ qA - f) J (6 . 5) 

It is interesting to observe that, in the limit ( 16.2b . our equation of motion ( 16.51 ) still contains terms describing "virtual" 
transitions to intrinsically excited states \aq) + \a$q). Since the two terms involving these transitions are of similar 
structure but of different sign, we do not expect a large contribution from them. This is no longer expected whenever 
the temperature T(q, t) is no longer small compared to the elementary excitations. In this case, the equation of motion 
(13.77b must be solved as it stands, together with the eq. ( 14.24b which determines the temperature T(q, t). 

One notices from Eq. (16.3b that the Boltzmann factor F a (q, f) becomes very small whenever the intrinsic excitation 
energy &' a (q) - £>' ao (q) is much larger than the temperature T(q, t). This means that, in the equation of motion ( 13.77b . 
only those intrinsic excitations \aq) need to be taken into account which correspond to non-negligible Boltzmann 
factors F a (q, t). 

The actual computing work concerns the calculation of the matrix-elements < H aa (q,q') and < H cp i aa ,(q,q''), on the 
one hand, and the numerical solution of the integro-differential equation (13.77b . on the other. One could convert the 
integral part of the equation into partial differential terms of the Fokker-Planck type by a moment expansion of the 
integral kernels. We do not think that this is useful, because the moment expansion introduces approximations and 
the resulting partial differential equation for R(q\,qi;i) is also difficult to solve numerically. Thus we think that a 
straightforward numerical solution of Eq. (13.77b on the basis of a discretization of the g.c. q is the best way to carry 
out practical applications of the theory. 
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Let us add that the techniques used in earlier calculations where intrinsic excitations were omitted Jgl |T(J could 
be employed also in actual applications of the present theory. For instance, in fission problems, it will probably be 
convenient to define a "scission region" beyond which fission fragments no longer exchange nucleons or internal 
energy and are subject only to Coulomb repulsion and decay. Then fragments characteristics such as mass-, charge- 
or kinetic energy distributions (prior to neutron and y-ray emission) could be evaluated from the values taken at large 
times by the solution of the equation of motion ( 13.77b obtained in the scission region or slightly beyond. 

7. Summary and Discussion 

a) Short review of what has been done 

The basis of our work is the method proposed by Hill and Wheeler of expanding the nuclear many-body state 
in terms of a set of basis functions which depend parametrically on deformation variables q (see Eq. (11.4b ) which are 
integrated over. In this way, the description of the nuclear shape and of its dynamical change as a function of time 
does not necessitate the introduction of superfluous degrees of freedom. 

This method which has been successfully applied to the description of single nuclear states has been applied to 
the nuclear density operator and its time evolution in the present paper. 

We used as basis functions states of independently moving quasi-particles defined as eigenstates <f> a (q) of the 
one-body part of H[{q) of the constrained Hartree-Fock-Bogoliubov (HFB) Hamiltonian (see Eq. ( 11. 101 )). It is a 
most welcome property of these states that, in the limit of vanishing collective large scale motion, they represent a 
self-consistently determined state of equilibrium at each given shape of the nucleus. 

A complication is produced by the non-orthogonality of these states. We coped with it by either using a biorthog- 
onal set of basis functions (see Eq. ( 11.15b and Section 13 or by introducing an orthogonalized version of the Hill- 
Wheeler states (see Eqs. (fSTTTh . dXT2l >). 

The basic assumptions of our theory are the ansatz ( 13.7b for the slowly time-dependent part pp(t) of the statistical 
operator and the Markov approximation (see text after Eq. ( 13.641 )) simplifying the equation of motion ( 13.641 ) for pp(t). 
The physical idea underlying the choice of pp(t) is that this slowly time-dependent part of the statistical operator 
depends neither on the detailed phase relations of the rapidly time-depending part of p(t) nor on the exact occupation 
amplitudes of the intrinsic excitations at a given shape of the nucleus which thus may be represented by a canonical 
distribution (see Eq. ( 13.81 )). 

In Section |U the temperature T which occurs in the canonical distribution is chosen by the requirement that the 
reference density pp{t) should correctly describe the total mean energy (see Eq. ( 14.1b ) and the energy at a given shape 
q of the nucleus (see Eq. ( 14.3b ) at all times. 

An important question is whether our theory implies a relation between the quantities describing the dissipation 
and fluctuation in the equation of motion (13.76b for pp(t). We dealt with this question in Section[5] 

b) Short discussion of what should still be done 

Obviously, a pertinent question is how one is to describe the emission of nucleons and photons during the fission 
process. We note that the emission during the dynamical process of fission is more difficult to treat than the emission 
from the fragments which have reached a thermal equilibrium. We are working on a description of emission processes 
during fission in the framework of our theory. The work will be submited for publication at our earliest convenience. 

A very important endeavour will be to perform actual numerical calculations on the basis of the presented theory. 
Indeed, it is one of the purposes of the microscopic theory to find out whether we are going to find a satisfactory 
agreement between calculated and measured results on the fission process. This is by no means obvious because not 
only the assumptions of our theory are at stake but also the very concept of whether it is at all possible to describe all 
the low energy phenomena of nuclear physics on the basis of one given effective nucleon-nucleon interaction. This 
implication is a much farther going assumption than the hypothesis of the existence of an energy-density functional 
enabling us to calculate ground state properties. 

It should be noted that calculations on the fission process at very low excitation energy have already been per- 
formed within the generator coordinate method without inclusion of intrinsic excitations [ 10] where an encouragingly 
good agreement between calculated and measured results has been found. The present theory would give the possibil- 
ity to extend these calculations to higher excitation energy and compare them with the existing careful experimental 
work. 
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Finally, a problem of its own merit would be to investigate in detail the classical limit of the theory presented. 
It surely is related to the ^-function limit of the narrow gaussian overlap of many-body wavefunctions pertaining to 
different values of the generator coordinates. This fact suggests that the Fokker-Planck approximation of our equation 
of motion ( 13.771 ) for the reduced density R{q\,q 2 ; t) might be the appropriate starting point for this investigation. 
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Appendices 

A. Some properties of the projection operators 

In Eq. (13.16b . we have introduced the matrix-representation of the projection P. As the temperature T depends on 
time, the operator P is unequal from and its matrix-representation is 

P a^a^qia-^q^a^qi, = & a\ClTp q^q^ q^qfi a^a^F ai {,C{2-> T{q 2 -, 0) (A.l) 

The matrix-representation of PP is thus obtained as 

a 3 q } a 4 q 4 

6 a i a 2 Sq i q 5 &qi ?6 "a$a 6 F a 2 

(q 2 , T(q 2 , f)) J] F b (q 2 , T(q 2 , t)) = (A.2) 

b 

where the r.h.s. vanishes because we have 

Y J F b (q 2 ,T(q 2 ,t))= 1 

b 

On the other hand, for the matrix-representation of PP, one finds 

F aiqiazqia^q-iatquP a 3 q 3 cnq 4 asqsd(,q(, — 
— & a\a 2 & q\q$& q^q^ asa^F ai 

(q 2 , T(q 2 , 0) ^ F h (q 2 , T(q 2 , f)) 

b 

= P a l q l a 1 q 1 a 5 q 5 a 6 q 6 (A. 3) 

The results 1A.2K and ( IA.3I ) can be written in the form 

P(i)P(t) = (A.4) 
P(t)P(t) = h) (A.5) 

respectively. As we have 



Q = -P (A.6) 

as a consequence of Eq. ( 13.171 ). we find 

P(i)Q(t) = = (l-Q(tj)&(i) (A.7) 

&(t)P(.t) = Q(t) = &(t)(l-Q(tj) (A.8) 
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from ( 1A.4I I. and ( IA.5t . Eqs. ( IA.7I ) and (IA.81 > imply the relations 



6(06(0 = 6(0 (A.9) 

£(02(0 = (A. 10) 

Q(t)P(f) = P(f) (A. 11) 

£(02(0 = (A. 12) 



For any given operator A, the matrix-representation of PA reads 

(aiqi\P(t)A\a 2 q2) = 2u P a x q i a 1 q 1 a i q 3 a i q l {a^\A\ a ^) 



a } q y a 4 q 4 

S ai a 2 F a2 (q2, T(q 2 , t)) ^{bq\\A\bq 2 ) 

b 

= S aiai A rea (q u q 2 ) — T(q 2 ,t) (A.13) 

where the "reduced part" A red {q\,q 2 ) of the matrix {a\q\\A\a 2 q 2 ) is defined by 

A red {q x ,q 2 ) := £<2>i<Zi|A> 2 <7 2 >«5 6lfr2 (A. 14) 

The operator P (t)A is thus seen to be diagonal with respect to the quantum numbers of intrinsic excitations and 
slowly time-dependent as the temperature T(q 2 , t) depends slowly on time. The factor dF a2 (q 2 , T{q 2 , t))/dT can be 
reformulated as in Eq. ( 13.46) and differs from a canonical distribution. 

B. Detailed derivation of the Green function ( 13.40) 

In what follows we show in detail how one obtains the solution ( 13.401 ) of the equation (13.32) . 

+ l -Q{t)tG{t, s) = Q(t)6(t - s) (3.32) 

This equation implies that the Green operator G(t, s) can be written in the form 

&{t,s) = Q{i)£{t,s) (B.l) 

Sunstituting dBTl into Eq. ( 13321 

&f)k<t, ff) + (3(0^-^ + l -rQ(t)LQ{t)k{t, s) = Q{t)5(t - s) (B.2) 

at n 

and using the relations (see Appendix lAl 

tit) = -P{t) = -Q(t)P(t), 
we find that the Eq. ( 1B.21 i can be written in the form 



2(0 



dK(t, s) x . i t A - 

P(f)^(f, s) + -LQ(t)K(t, s) - 6{t - s) 



dt ' y7 ' h 

The operator ^T(f, s) is thus seen to satisfy the equation 

dK{t, s) 



(B.3) 



- W(t)K(t, s) = <5(f - s) (B.4) 



5? 
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with 

W(t) := P(t) - l -LQ{t) (B.5) 

n 

For solving (IB.4t in the range f > s, we put 

£(f, s) = o (t - s)K(t, s) (B.6) 

with the initial condition 

K (t, t) = 1 (B.7) 

Substituting ( IB. 6b into ( IB ,4b yields 

- s)k(f, s) + e (t - s) - *)| = - (B.8) 

and, because of 

(J(f - s)K(t, s) = 5(f - s)l(f, f) ^ 5(f - s), (B.9) 



we find the homogeneous differential equation 

df 0, s) 



= W(f)£(f,s) (B.10) 



5? 

The equation ( IB. 10b and the initial condition ( IB. 71 ) show that A^f, s) can be written in the form 

K(t, s) = 1 + J dr W(t)K(t, s) (B.l 1) 
The integral equation ( IB .lib has the solution 

k(t,s) = tL£ dT * ( - T A (B.12) 



where f is the time-ordering operator. The form JB. 12b is thus seen to be a short way to formulate the iterative solution 
of Eq. dBTfl , Finally, from (ED, dB~6b . and dB~T2l we find the form 

6{t, s) = 8 (t - s)Q{t)t U, dT (B.13) 

of the Green operator. 

C. Comments on the non-hermitian nature of pp(t) 

As we have already pointed out in Section [3] our ansatz ( 13.7b for the approximate form of the density operator is 
not hermitian. The reader might ask the justified question whether it would not be preferable to derive and solve an 
equation of motion for the hermitian ansatz 

p H P (t):=i(pp(t)+pl(t)) (C.l) 

rather than solving the equation of motion ( 13.76b for p pit) and subsequently use its hermitian part for the evaluation 
of physical quantities. 

If we use the orthogonal basis functions v F fl (^) (see Section|2|, the matrix representation of p h p (t) has the form 

QVaM^P^aMi)) = ^(F ai (quT(q l ,t)) + F a2 (q 2 ,T(q 2 ,t)j)R(q u q 2 -,t) (C.2) 
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It is seen that ( IC 2b is close to the matrix representation ( 13.71 ) of pp(t), if the Boltzmann factors F ai (q\, T{q\, t)) and 
F a ^(q 2 , T(q2, f)) do not differ very much one from another for pairs of values of the g.c. where the reduced density 
matrix R{q\,q 2 ; t) is substantially different from zero. We expect that the cases where (IC.21 > differs substantially from 
the non-hermitian case 

0F ai (qi)\pp(t)\*¥ a2 (q 2 )) = 8 aia2 F 02 (q 2 , T(q 2 , t))R{q u q 2 ; t) (C.3) 

are more important for small excitation energies (and temperatures), where the canonical ansatz is anyhow less well- 
founded. 
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